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Предложены быстродействующие алгоритмы деления полиномов в арифметике по модулю два, позволяющие работать сданны- 
ми в параллельном виде: модификация алгоритма одновременного определения старших и младших разрядов частного, а также 
разработанный авторами матричный алгоритм деления полиномов. Приведён пример реализации на ПЛИС фирмы АІІега кодека 
двоичного циклического помехоустойчивого кода, использующего матричный алгоритм ускоренного деления полиномов. 


Использование помехоустойчивых кодов счита- 
ется основным средством обеспечения требуемой 
достоверности передачи данных. В частности, 
двоичные циклические помехоустойчивые коды 
используются в системах промышленной автома- 
тизации, радиоинтерфейсах, кодеках видеопосле- 
довательностей, шинах передачи данных процессо- 
ров. В связи с увеличением скорости передачи дан- 
ных требуются быстродействующие помехоустой- 
чивые кодеки, которые смогут работать на частоте 
шины передачи данных. 

Основной операцией при кодировании и деко- 
дировании двоичным циклическим помехоустойчи- 
вым кодом является операция деления полиномов в 
арифметике по модулю два, поэтому необходимы 
быстродействующие алгоритмы деления. Обычно 
полиномы представляются следующим образом: 

п - 1 т-1 

А (х) = (2Х-Ы*" -1- '); в(х) = С^ь т _ н х "" н ), (і) 

і=0 7=0 


где а„ й,е{0,1}; /=0,я-1;у=0,т— 1; гі>т; п - количе- 
ство разрядов делимого, т - делителя. 

В [ 1] предложен алгоритм деления полиномов в 
арифметике по модулю два и доказана теорема для 
математического обоснования основной идеи ал- 
горитма (одновременное определение старших и 
младших разрядов частного), формулировка кото- 
рой представлена ниже. 

Теорема 1 . Если полиномы вида (1) делятся без 
остатка, то полиномы, полученные из исходных 
путем изменения порядка следования коэффици- 
ентов на обратный (А' (х)=х"А(хг 1 ); Е(х)=зСВ{х~ 1 )), 
также делятся без остатка. Полученное таким обра- 
зом частное будет иметь обратный порядок следо- 
вания коэффициентов по сравнению с частным от 
деления исходных полиномов (С(х)=х г С(х^')). 

В данном алгоритме автором было наложено 
ограничение: младшие разряды делимого и делите- 
ля должны равняться 1. Для устранения ограниче- 

19 


Известия Томского политехнического университета. 2006. Т. 309. № 1 


ний алгоритма нами было сформулировано и дока- 
зано следствие из Теоремы 1 [2]. 

Следствие из Теоремы 1. Полиномы вида (1) де- 
лятся с остатком, если п-р<т, где р — количество 
младших незначащих нулей делимого. 

Доказательство (от противного). 

Обозначим через р количество младших незна- 
чащих нулей полинома А. Тогда, 

А(х) = +... + а р х р ) = 

= (Ъ^хГ- 1 + ...+ Ь 0 х 0 )(С„_ т х"- + ... + С р х р ) = 

= 5(х)С(х); 

+... + а 0 х°) = 

= (Ъ т _ 1 х т -'+... + Ъ 0 х°)(С„_ т х р - + ... + С р х р ); 
(а п _,_ р х"- 1 - р +... + а 0 х°) = 

= (ъ т _у'~' +...+ъу\с, 1 _ р _у- р - т +...+с 0 А 

если т>п-р, то степень п-р-т<0. 

Что и требовалось доказать. 

Примечание: предполагается, что младший раз- 
ряд делителя ненулевой, иначе степень делителя 
вычисляется без младших незначащих нулей. 

На основании Теоремы 1 и Следствия из Теоре- 
мы 1 разработан «Двусторонний алгоритм ускорен- 
ного деления полиномов в арифметике по модулю 
два», при котором одновременно определяются 
старшие и младшие разряды частного [2]. В отли- 
чие от алгоритма [1], данный алгоритм является 
универсальным (в нём не накладывается никаких 
ограничений на делимое и делитель) и всегда обес- 
печивает корректное вычисление частного. 

Пусть п - разрядность делимого, а т - разряд- 
ность делителя, тогда разрядность частного равна 
г=т-п+ 1 ; А — делимое, В - делитель, С - частное, 
а Б - остаток. 

Двусторонний алгоритм ускоренного деления по- 
линомов в арифметике по модулю два. 

Начало: даны полиномы вида (1). 

Шаг 1: Подсчитываем количество младших нез- 
начащих нулей делителя, обозначаем его через р. За- 
поминаем р последних разрядов делимого и отбрасы- 
ваем их. Получаем новое делимое А { =а п _ х ...а=а\ л ...с^. 

Частное С І =С 5 ...С^, где х=п-т, С,= 0, /=0,і. По- 
лагаем номер шага /'= 1 , г,=0, /,=$. 

Шаг 2: Если Д^О, подсчитываем количество 
старших незначащих нулей Ь, и количество млад- 
ших незначащих нулей К, в делимом А=а[ і _ ] ...а 1 !І . 
Полагаем //=/- Ц; /•/=/} +Д-. Отбрасывая в Д старшие 
и младшие незначащие нули, получаем новое дели- 
мое Д разрядности Иначе идём на шаг 5. 

Шаг 3: Если число разрядов п>2(т—р), форми- 
руем двоичное число (новый делитель): 

В =Ь , ...Ъ 0...0& . ...Ъ , 

і т-1 р і_ ѵ _1 т-1 р ’ 


разрядности Д и Д совпадали). Получаем двоичное 
число Д+1=ДѲД. Заносим в разряды Ди //частного 
единицы (полагаем С, ; =С с =1, /,+1=//-1; г ( +1=ДН; 
/=/+ 1). Идём на шаг 2. Иначе идём на шаг 4. 

Шаг 4: Формируем новое делимое А следую- 
щим образом: приписываем в конец Д(гД1) нулей 
и р разрядов, запомненных на шаге 1 : 

Д =а і „_ 1 ... а‘ р Ф..0 а р ...а 0 . 

І} + 1 

Делим А на В обычным способом. Обозначим 
через /частное отделения, О - остаток / =В •/+/). 

Частное от деления / на В формируем следую- 
щим образом: С=СѲ/ 

Остаток от деления равен Г). Идём на шаг 6. 

Шаг 5: Частное от деления / на В равно С. Ос- 
таток отделения і)=а г ..й 0 . Если остаток 0=0. то по- 
линомы делятся нацело. 

Шаг 6: Конец. 

Исследовав особенности построения двоичных 
циклических помехоустойчивых кодов, заметив 
полезные свойства при построении кода с помо- 
щью образующей матрицы в систематической фор- 
ме, мы разработали «Матричный алгоритм уско- 
ренного деления полиномов в арифметике по мо- 
дулю два» [3]. 

Введём обозначения 

А(х ) = А'(х) ® А(х); А'(х) = л _ Ьі х“ _1 ”‘ , 

і=0 

^ л-1 т— 2 

А(х) = ^ а іі і і х п ~'~ і и В(х) = ^. 

і=п-т+1 ^=0 


Определим матрицу /: 



, где ^(х)=К вм (х" 1 ') (остаток от деления 


2 


п- 


т 


хГ'~‘ на делитель В{х)), іе (0,п-т). 

Матрицу У определим следующим образом: 



, где у,(х) 


В(х) 


(частное от деления 


[_Уп- 


т 


х"“ 14 на делитель В(х)). 

Матричный алгоритм ускоренного деления поли- 
номов в арифметике по модулю два. 

Начало: Даны полиномы вида (1), где А - дели- 
мое, а В - делитель. 

Шаг 1: Вычислить матрицу /для полинома В(х): 


/ = 


2 0,т-2 

г 1,т-2 


2 


п-т,т-2 


где ц, - число нулей, которые определяются из уело- Шаг 2: вычислить матрицу У для делимых раз- 
вия: т-р+ 1 +д,+т-р + 1 =п, (число нулей таково, чтобы рядности п: 
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У = 


У О.п-т 

Уі,п-т 

У п—тм—т 


Уо, о 
Уі о 

У п-т ,0 


Шаг 3: Вычислить остаток А ' • 2® А . 


Шаг 4: Вычислить частное А' • У. 

Шаг 5: Конец. 

Для математического обоснования верности ал- 
горитма сформулируем и докажем ниже предста- 
вленную теорему [3]. 

Теорема 2. Если А' , 2®А = 0, то А(х) делиться на 
В(х) без остатка, иначе остатком будет являться ре- 
зультат. 

Доказательство. 

Рассмотрим разложение х" -1- '' на В(х): 

х п ~ 1 ~‘ =В{х)-у і ® 2 і {х), (2) 

где і=0,п-т\ і, - і - ая строка матрицы I. 
Выражение (2) подставим в А'(х ): 

п-т 

А '(х)=Х а «- 1-Г •к"''"' = 

1=0 


= [X а «-\-і ’УіУ В ( х ) ф X а «-и • 2 , С*) ; 


А(х) = А’(х) 0 А(х) = ->Дх 


<В(х)®[^а п _ ) _ Г 2 і (х)® X ■х п - 1 ~ і 1 (3) 


Из (3) видно, что остаток равен 

2а-і-г 2 ,- (*)© X 


что соответствует 

и-т й-1 

И'.2®й=Х«„-ы- 2 ,'® I <* п -ыХ п -'- 1 - 

/—О і -н -т + 1 

Что и требовалось доказать. 

Следствие из теоремы: 

И' • У = ^ а п -і-іУі , А'У= частное. 

/=0 

Для быстрого формирования матриц 2 и У, 
предложим и докажем следующие выражения: 

1. 2 ,(х) = 2 і+І (х)-х® 2 /+І/п _ 2 -В(х), где _ 

старший разряд с,+ь і=п—т—\ ,0; 

2 - Е,(-Д = К + і(-Д--г® 2 . + і / „_ 2 , где /=й— м-1,0. 

Доказательство методом математической ин- 
дукции: 

Базис. і=п-т. Тогда 

1 п- т (х)=В т (Г-')=В(х) и у п _= 1 , что очевидно. 
Индукция. ()</</. Тогда 


х у+1 =В{х)-у^{х)® 2 і+1 (х)\ 

X ‘ = В(х)- у ]+] (х)-х® 2 /+ , (х) ■ X = 

= В(і) ■ у і+1 (і) • і Ѳ 2 ;+1 (х) ■ і Ѳ ^ -В(х)®Ѵ] ■ В(х) = 
= В(х)[у 1+] (х) ■ х 0 ѵ,. ] © [г у . +1 (х) • х® у, • Я(*)1- 
Ч = 1 еС ЛИ ^■+1, т _2 = 1=>Т = ^ +1т _2. 

Отсюда следует, что у,(х)=у ;+1 (х)-хѲг ;+ і, и _ 2 и 
^)=^ +1 (х)-х®г у+ і, т _2-5(л:). 

Что и требовалось доказать. 

Для определения быстродействия алгоритмов 
деления полиномов в арифметике по модулю два 
поставлен компьютерный эксперимент, в котором 
использовались: 

• персональный компьютер: Іпіеі Сеіегоп 1 ГГц, Іп- 
іеі 815ЕРВ сЬірзеІ, ВОКАМ 512 Мб, ИБО 60 Гб; 

• операционная система Місгозой \Ѵіпс1оѵѵх ХР 8Р2; 

• разработанный программный эксперименталь- 
ный комплекс [4]. 

Эксперимент проводился над 30-разрядными 
кодовыми словами. Использовался образующий 
полином вида: х 6 +х 4 +х 3 +х 2 +х 1 + 1 . 

Результаты компьютерного эксперимента све- 
дены в таблицу. 


Таблица. Результаты компьютерного эксперимента 


Количество ко- 
довых слов, тыс. 

Время работы алгоритма, с і 

Стандартный 

Двусторонний 

Матричный 

і 

0,110 

0,110 

0,080 

5 

0,410 

0,371 

0,220 

10 

0,791 

0,761 

0,471 

15 

1,121 

1,062 

0,651 

20 

1,422 

1,372 

0,822 

25 

1,763 

1,752 

1,082 

30 

2,113 

2,003 

1,222 

35 

2,413 

2,363 

1,523 

40 

2,774 

2,704 

1,773 

45 

3,064 

2,944 

1,923 

50 

3,425 

3,235 

2,053 


Из таблицы видно, что двусторонний алгоритм 
в среднем на 4 % быстрее, чем стандартный алго- 
ритм [1], а матричный алгоритм быстрее стандарт- 
ного в среднем на 40 %. 

Преимущество двустороннего алгоритма про- 
является при большой длине кодового слова, когда 
не хватает памяти в системе для размещения ма- 
трицы остатков матричного алгоритма. В практи- 
ческих задачах рекомендуется применять быстро- 
действующий матричный алгоритм деления поли- 
номов в арифметике по модулю два. 

В качестве примера нами реализован кодек 
двоичного циклического помехоустойчивого кода с 
использованием матричного алгоритма деления на 
ПЛИС фирмы Мега семейства Асехі . Для исполне- 
ния на ПЛИС выбрана структура кодека [5], рис. 1. 
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На рис. 1 обозначены: 

• с\ѵ[п- 1 ..0] - входные данные (кодовое слово для 
декодирования или информационные символы 
для кодирования); 

• сЬ[г- 1 ..0] - контрольный блок; 

• Х[г-1..0] - синдром ошибки; 

• іЪ[к- 1..0] - информационный блок; 

• соіес - входной сигнал, который переключает 
кодек в режимы кодирования и декодирования; 

• ре : - позиция /-ой ошибки; 

• /„■ - ис прав ляющая комбинация для /-ой ошиб- 
ки, /е(1,се); 

• ст[п- 1..0] - выходные данные кодека, т.е. ко- 
довое слово, если сосІес= 1 и информационный 
блок с сигналом ошибки, если соіес= 0, где п - 


разрядность кодового слова, г - разрядность 
контрольного блока, к - разрядность информа- 
ционного блока. 

Логические блоки: 

• КС 1 - формирования контрольных разрядов с 
использованием быстродействующего матрич- 
ного алгоритма деления полиномов в арифме- 
тике по модулю два; 

• ПЗУ 2* - вычисления позиций ошибок; 

• КС 2* - формирования исправляющей комби- 
нации. 

В качестве двоичного циклического помехоу- 
стойчивого кода был выбран код БЧХ (31,21), ко- 
торый используется в стандарте формата радиовы- 
зова РОС8АС. Образующий полином данного кода 
равен 8(х)=х' 0 +х 9 +х*+х 6 +х 5 +х ! + 1 . 



Рис. 1 . Структура кодека двоичного циклического помехоустойчивого кода 
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Кодек двоичного циклического помехоустой- 
чивого кода спроектирован с помощью языка опи- 
сания аппаратуры Ѵегі1о§. Структура программы 
кодека представлена на рис. 2. 

сосІес_рос8ад«Ь © — >сазе_рос8адк :ту_са®е2 

»са8е_рос5адз :ту_саяе1 

* гот_рос5ад!5 :ту_готЗ 

> сііѵ_рос»адЬ :ту_сііѵ 

Рис. 2. Структура программы кодека 

Работа кодека промоделирована в пакете 
Мах+рІшІІ фирмы Акега. Результаты моделирова- 
ния для ПЛИС ЕРР10К30ЕТС144-1 представлены 
на диаграммах в шестнадцатеричной системе исчи- 
сления, рис. 3, 4. 

На диаграмме представлены следующие сигналы: 

• Даш - входные данные кодека (кодовое слово 
или информационный блок); 

• У 1х_с1а(а - исправленное кодовое слово; 

• геасіу — сигнал готовности кодека; 

• сіоск - синхронизация; 

• (Ііѵ\сІіѵісІепёН - информационный блок в кодо- 
вом слове; 


• ёіѵ\<ИѵісІепёЬ — контрольный блок в кодовом 
слове; 

• сосіес - сигнал, переключающий кодек между 
режимами кодирования и декодирования; 

• у_с!іѵ\гетіпсІег - контрольные разряды или син- 
дром ошибки (зависит от сигнала соДес); 

• _са$е 1\/іх_егг - исправляющая комбинация для 
первой ошибки; 

• _са$е2\рх_егг - исправляющая комбинация для 
второй ошибки; 

• еггог - сигнал об обнаружении неисправимой 
ошибки. 

Диаграмма моделирования работы кодека в ре- 
жиме исправления ошибок в кодовом слове пред- 
ставлена на рис. 4. Ошибки находились в двух стар- 
ших разрядах кодового слова. 

В заключение необходимо отметить, что разрабо- 
танный быстродействующий алгоритм матричного 
деления полиномов в арифметике по модулю два по- 
зволил создать быстродействующий кодек двоичного 
циклического помехоустойчивого кода. Время рабо- 
ты кодека для ПЛИС ЕРЕ10К30ЕТС144-1 в режимах: 

• кодирования 10 нс; 

• декодирования 15 нс. 


ІМате: 

ѴаІие: 

15.0(13 

20.0пз 25.0 пз 

ЗО.Опз 


Зб.Опз 
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00ЕАВ400 

-к# 1іх_с1а1а 

Н 00000000 


00000000 
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000 
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1 
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о 

о 
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о 
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о 
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00001000 
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о 
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о 
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о 

о 

X 
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Рис. 3. Диаграмма моделирования работы кодека в режиме кодирования 
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Рис. 4. Диаграмма моделирования работы кодека в режиме декодирования 
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